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We study the intrinsic friction of monolayers adsorbed on solid surfaces from a gas phase or vapor. 
Within the framework of the Langmuir model of delocalized adsorption, we calculate the resistance 
offered by the mobile adsorbate's particles to some impure tracer molecule, whose diffusive random 
motion is biased by a constant external force. We find that for sufficiently small driving forces 
the force exerted on the tracer shows viscous-like behavior. We derive then the analog of the 
Stokes formula for two-dimensional adsorbates, calculate the corresponding friction coefficient and 
determine the stationary particle distribution in the monolayer as seen from the driven impurity. 
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When an ambient gas phase or vapor is brought in con- 
tact with a solid, some portion of it becomes reversibly 
attached to the solid surface in form of an adsorbed layer. 
Such layers are important for various technological and 
material processing operations, including, for instance, 
coating, gluing or lubrication. 

Following the work of Langmuir, equilibrium proper- 
ties of the adsorbates have been extensively studied and a 
number of significant developments have been made |Q. 
As well, some approximate results have been obtained 
for both dynamics of an isolated adatom on a corrugated 
surface and collective diffusion, describing spreading of 
the macroscopic density fluctuations in interacting ad- 
sorbates being in contact with the vapor ^J,|| . 

Another important aspect of dynamical behavior con- 
cerns tracer diffusion in adsorbates, which is observed 
experimentally in STM or field ion measurements and 
provides a useful information about adsorbate's viscosity 
or intrinsic friction. In this regard, analysis of tracer dif- 
fusion is not only a challenging question in its own right, 
but is also crucial for understanding of various dynamical 
processes taking place on solid surfaces. To name but a 
few, we mention spreading of molecular films on solid sur- 
faces ||], spontaneous or forced dewetting of monolayers 
Q or island formation j6|. However, apart of a slightly 
artificial ID model , available studies of tracer diffu- 
sion in adsorbed layers focus on strictly two-dimensional 
(2D) situations (see, e.g. Refs. |8|~|i"o|), excluding the 
possibility of particles exchanges with the vapor. 

In this Letter we present first results on intrinsic fric- 
tion in 2D adsorbed monolayers composed of mobile par- 
ticles undergoing continuous exchanges with the vapor 
phase. The system we consider here corresponds to the 
generalized Langmuir model of adsorption with adsor- 
bate lateral diffusion |l| and consists of three key ingredi- 
ents: (i) a solid surface containing some concentration of 
adsorption sites, (ii) a vapor phase and (iii) a monolayer 
of adsorbed hard-core particles, which perform activated 



random motion between the adsorption sites subject to 
the hard-core exclusion interactions (one particle per ad- 
sorption site at most), and undergo continuous exchanges 
(desorption/adsorption) with the vapor phase (Fig.l). 




Fig. 1. Adsorbed monolayer in contact with a vapor. 
Grey spheres denote the monolayer (vapor) particles; 
the smaller black sphere stands for the driven tracer 
particle. 



To determine the intrinsic friction of such an adsorbate 
we wish to probe the resistance offered by the monolayer 
particles to some external perturbance. To this purpose, 
we add an extra hard-core particle - a tracer, which may 
move along the surface only and can not desorb. Next, we 
suppose that the tracer is subject to an external force F, 
which favors its motion in the preferential direction such 
that it ultimately attains a constant velocity V. Our aim 
is then to evaluate the velocity- force relation V{F) in the 
limit of small F and to determine the stationary particle 
distribution in the perturbed monolayer. General results 
for arbitrary F will be presented elsewhere. 
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More precisely, our model is defined as follows (Fig.l): 
We suppose first that the adsorption sites form a regular 
2D square lattice of spacing a. Occupation of the lat- 
tice sites is characterized by the set of time-dependent 
variables rjft, f being the lattice vector of the site with 
Cartesian coordinates (x,y); rjp equals 1 if the site r is 
occupied by an adsorbed particle and zero otherwise. 

Next, we suppose that the particles from the vapor, 
(considered as a reservoir, which is maintained at a con- 
stant pressure), may adsorb onto any vacant lattice site 
at a fixed rate f/r*. Then, the adsorbed particles may 
move randomly by hopping with a rate I /At* to any 
of four neighboring lattice sites, which process is con- 
strained by hard-core exclusion preventing multiple oc- 
cupancy of any site, or desorb from the lattice at a site- 
and environment-independent rate g/r*. For simplicity 
of exposition, we assumed here that typical adsorption, 
desorption and jump times are equal to each other; these 
times, i.e. T a d, Td es and Tj ump , respectively, can be read- 
ily restored in our final results by the mere replacement 
/ -* T*f/r ad) g -> T*g/r des and I -> It* /t^ 

The tracer particle is initially put at the lattice origin. 
We stipulate that this only particle is not allowed to des- 
orb and in addition to random thermal forces experiences 
an action of a uniform external force F = (F,0). Dynam- 
ics of the tracer is then defined via standard rules [ pT[ : 
the tracer, which is at the site R, R = (X, Y), at time 
t, waits an exponentially distributed time with mean r, 
(which is in the general case different from r*), and then 
attempts to hop onto one of the neighboring sites, R', 
with probability p(R\R'), defined as 



p(R\R>) = exp (t{Fe v )) / £ exp(^{F^)) , (1) 



where e v = R' — R, v,p = l,2,..,z, z being the coor- 
dination number of the lattice, while [3 is the reciprocal 
temperature. The hop is fulfilled only if the target site 
appears to be vacant; otherwise, the tracer remains at R. 

Now, several comments are in order. First, we note 
that particles diffusion does not matter in case when the 
tracer and the inter-particle attractive interactions are 
absent; here, the monolayer is homogeneous and the local 
occupation variables tend as t — > oo to the same value 
Pl = 1/(1 + g/f), the latter relation being the Langmuir 
adsorption isotherm [jD. In the presence of an impure 
molecule, whose hopping rates differ from that of the 
monolayer particles, particles diffusion is significant and 
couples effectively the occupation variables of different 
lattice sites. This results, as we proceed to show, in the 
appearence of rather complex density profiles. Second, 
we remark that in the model under study the number 
of adsorbed particles is not explicitly conserved. The 
conserved particles number (CPN) limit can be achieved 
in the stationary regime by setting / and g to zero, while 
keeping their ratio fixed, f/g = Pl/0- ~ Pl)- This limit 
corresponds to standard models of tracer diffusion in 2D 



hard-core lattice gases. Lastly, we remark that our model 
can be thought off as a certain generalization of a model 
for dynamic percolation proposed in ]l2] . 

Let P(R, rj) denote the probability of finding at time t 
the tracer at the site R and all adsorbed particles in the 
configuration r/ — {r/p}. Further more, let rf ,v denote 
the configuration obtained from r\ by the Kawasaki-type 
exchange of the occupation variables of the two neigh- 
boring sites r and r + e v , and ff be the configuration 
obtained from rj by the replacement rjp — > 1 — rjp, which 
corresponds to the Glauber-type flip of the occupation 
variable due to the adsorption/desorption events. Then, 
counting up all events which can result in the configu- 
ration (R,rj) at time t or modify it, we write down the 
following master equation 

d t P(R, V ) = [p(0\e v )P(R-e v ,r,)(l- m ) - 

-p(0\e u )P(R,r,)(l- m+s j] [P(R,V P 'l- 

- P(R, v )] + [(i - v?)P(R, f) - vfP(R, v) 



+ 



[vfP(R,f)-(i-Vf)P(R,v) 



(2) 



where the symbol r under the summation sign means 
that the sum runs over all the lattice sites. 

Now, we are in position to obtain the X-component of 
the tracer velocity V. Multiplying both sides of Eq.(||) 
by X and summing over all (i?, r/), we have 



a 

V x = - 

T 



p(o|M)) -p(o|(-M)) x 



X 



l-fc(-a,0)jj, k(X) = J^ m+x ,P(R, V ) (3) 



where k(X) is the probability of having at time t an ad- 
sorbed particle at position A = (X x ,X y ), defined in the 
frame of reference moving with the tracer particle. In 
other words, k(X) can be interpreted as the density pro- 
file as seen from the moving tracer. 

Consequently, Vx depends on the particle density in 
the immediate vicinity of the tracer. Note, that if the 
monolayer is perfectly stirred, i.e. k(X) = p^, (which 
implies decoupling of the occupation variables), we would 
obtain from Eq.(S) a trivial mean- field result Vx = Vq = 
a(p — q)(l — pl)/t, which states that the tracer jump 
time r gets merely renormalized by a factor (1 — p^) -1 , 
which defines the frequency of successful jump events. 
However, this is not the case and A;(A) ^ pl, except for 
|A| — > oo. Moreover, k(X) is a complicated function of 
Vx, which renders Eq.(0) to be strongly non-linear. 

Equations describing the time evolution of fc(A) can be 
found from (||) by multiplying both sides of it by il^ +X 
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and summing over all configurations (i?, 77). In doing so, 
we find that these equations are not closed with respect 
to fc(A), but are coupled to the third-order correlations, 



T(A, e u ) = Vr + xVr +Sv P(R, V) 



R,7) 



In turn, if we proceed further to the third-order corre- 
lations, we find that these are coupled respectively to 
the fourth-order correlations. Consequently, in order to 
compute Vx , one faces the problem of solving an infinite 
hierarchy of coupled equations for the correlation func- 
tions. Here we resort to the simplest non-trivial closure 
of the hierarchy in terms of k{X) representing T(A, e„) as 



(5) 



We hasten to remark that the decoupling in Eq.(||) 
provides exact results for the analogous ID model in the 
CPN limit |l3Ur^ 1 and serves as a very good approxima- 
tion for the ID model with a reservoir |?]]. We set out to 
show in what follows that in the CPN limit our results 
reproduce exactly the classic results of Nakazato and Ki- 
tahara P], which are exact for pL <C 1 or Pl ~ 1, and 
serve as a very good approximation for any intermediate 
Pl (§]■ Since adsorption/desorption processes are essen- 
tially linear, we expect that such a closure will provide 
an accurate description for arbitrary / and g. 

Employing the approximation in Eq.(|5|), we obtain for 
k(X) the following closed evolution equations 



d t k(X) = L fc(A) 
f + 9 , 1 



4t* 



+ -$^p(0|£,)(l -*(?„)) V*„, 



(6) 



Equation (g) holds for any site A, except for the four 
sites in the immediate vicinity of the tracer: k(X) at these 
four neighboring sites is essentially perturbed due to the 
asymmetric hopping rules of the tracer and obeys 



k{S v ) = (l fc(A)) 



+ -p(Q\e v )k{2e u ){l-k(e„) 

A=e„ T 



Solution of Eqs.(||) to (||) in the limit t — > 00 can be 
readily obtained by applying the generating function ap- 
proach and reads: 



(4) fe(A) = p L + 



2-KjOtp 



-Pl 4,o - A. 



a A) 



(9) 



where Ag v = l+±T*p( Q\e v ) (l-k( e v )) /It, a F = J2 V ^e„ + 
Ag/l(l- pl), = 2(^/4,0^,0 + 4^)^, while 



A 



and 



4 



4 



Aj,/2<7 -7T 



d 2 k- 



^exp 



i{kX) 



v /4,o-4-<r,o cos(ak x ) + A 0>(7 cos(ak y ) 

1/ 4,0^4-CT.O + 4,<T 



(10) 



(11) 



Lastly, we have to determine four particular values k(e v ) 
which appear on the rhs of Eq. (^|) . This can be done by 
setting A = e v , v = 1,2, ..,z, to the lhs of Eq.(^|) and 
solving the resulting system of four linear equations. In 
doing so, we define k(e„) explicitly, as functions of Vx 
and of the characteristic microscopic parameters. 

Consider now the behavior of the tracer velocity Vx in 
the limit aj3F <C 1. First, we have from Eq.([|) that 



V x ~ — {PaF(l 



Pl) - 5k} . 



(12) 



where 5k = k(a,0) — a, 0) denotes the density jump 
in the vicinity of the tracer. Note, that since 5k > 0, the 
tracer velocity Vx is always smaller than the mean-field 
prediction Vq. Next, from Eqs.(0) to (O) we find 



5k 
aj3F 



2 PL (i-p L y 3 P 



1 It WA afl {F = Q) 



Ak{X) = Vv e ^fc(A) = V \k(X + e v )-k(X)\ (7) 

u Jo W 



exp 



a Q w 



24,o(F - 0) 



h{w) I (w) 



(13) 



I n {w) being the modified Bcsscl function of order n. 

Consequently, we find from Eq.(|l2|) that in the limit 
a(3F <C 1 the force-velocity relation attains the physi- 
cally revealing form of the Stokes formula, i.e. F = (Vx, 
which signifies that the friction force exerted by the 
monolayer particles on the driven tracer is viscous. The 
corresponding friction coefficient £ is given by 



- -k{e v )p{e v \Q)(l - k(-e u )) - -Lv_^fc(0) + 4 (8) 
r V / 4r* t* 

Note that Eqs(||) to (g) constitute a closed system of 
equations, which suffice in principle the computation of 
the density profiles and of the tracer velocity. However, 
these equations are non-linear, since k(e v ) enters the pref- 
actor before the gradient terms, which makes such a com- 
putation to be a non-trivial problem. Below we consider 
the stationary solution of Eqs.(||) and (||). 



C Co ~t~ Ccoopi 

where Co = 4t/<t 2 /3(1 — pl) and Ccoop follows 



c 



coop 



Co 



(5k/al3F) 



I- PL- (5k/a(3F) ' 



(14) 



(15) 



(5k/a(3F) being defined in Eq.(|lJ). While Co is the 
typical mean-field-type result (see the discussion follow- 
ing Eq.(0)), the second term has a more complicated 
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origin and is associated with the cooperative behavior: 
de-homogenization of the monolayer by driven impure 
molecule and formation of stationary density profiles, 
whose characteristic parameters depend on the velocity 
Vx- This second term can be small when either pl <C 1, 
or r* — > (perfect mixing), but dominates the overall 
friction for any intermediate values of systems parame- 
ters. We finally remark that in the CPN limit Eq.jl^) 
reduces to the classical result obtained in ||, which is 
known to be exact when pi, <C 1 or ~ 1 and rep- 
resents a very good approximation for any intermediate 
value oi pl || • 

Lastly, we discuss the characteristic features of the 
monolayer density profiles as seen from the stationary 
moving tracer, Eq.(|9|), at large distances in front of and 
past the tracer. It follows from Eq.(j^) that in the limit 
| As, | — > oo, (X y = 0) the density profiles obey 



k(X x , 0) ~ p L ± y (^) 1/Z exp( ± \^ln(z ± )) , (16) 

where the sign "+" ("-") corresponds to the domain X x > 
(X x < 0), while z± are two of four eigenvalues of the 
operator L in Eq. (|6|) , which are given by 



1/2 



A, 



.1 (OLF . 



i±wi- 



2A a oA- a o 



(17) 



Hence, in the domain A^ > the monolayer density is 
higher than pl, which means that there is a "traffic jam" 
type region in front of the tracer, which impedes its mo- 
tion. This region decays exponentially with the distance. 
In turn, past the tracer there is an exponentially decay- 
ing depleted region in which the density is lower than pl; 
since Z- < z+, the depleted region is more extended in 
space than the condensed one in front of the tracer, such 
that the density profiles are asymmetric with respect to 
X x = 0. It is interesting to note that for k(X) as in Eq.(j|) 

the sum yiy(fcfA) — Pl) = 0, which means that the tracer 
does not perturb the global balance between adsorption 
and desorption. This is not, however, an a priori evident 
result in view of the asymmetry of the density profiles. 

The salient feature of the behavior past the tracer is 
that in the CPN limit Z- = 1 at any fixed F, which 
signals that Eq.([l^) is no longer valid and correlations 
get somewhat stronger. Indeed, we find that 



k(X x ,0) - pl ~ -const x | X x | 3/2 , X x 



-oo 



(18) 



i.e. the correlations fall off with the distance as a 
power — law\ Remarkably, this implies that in the CPN 
limit the mixing of the monolayer is not at all efficient and 
there are considerable memory effects - the host medium 
remembers the passage of the tracer on a large space and 
time scales. This situation can be realized experimentally 
for ultrathin liquid films confined in narrow slits between 
solid surfaces, e.g., in boundary lubrication. We note 



also that the power-law behavior in Eq.([L8|) can be ob- 
served as an intermediate scale decay for adsorbed layers 
exposed to a low vapor pressure. 

To conclude, we have studied analytically the intrin- 
sic frictional properties of 2D adsorbed monolayers, com- 
posed of mobile hard-core particles undergoing continu- 
ous exchanged with the vapor. By analysing the force- 
velocity relation of a driven impure molecule - a tracer, 
which is designed to probe the resistance offered by ad- 
sorbate particles to the external perturbance, we have 
derived the analog of the Stokes formula for 2D mobile 
adsorbates and calculated explicitly the corresponding 
friction coefficient. Besides, we have determined the sta- 
tionary density profiles, which emerge in the adsorbate 
in response to the presence of a driven impurity, as well 
as obtained explicit results for both the density jump in 
the immediate vicinity of the tracer and the asymptotical 
density relaxation at large separations from it. 
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